3site Multisubstrate-Bound State of Cytochrome P450cam.

We identified a multisubstrate-bound state, hereby referred as a 3site state, in cytochrome P450cam via integrating molecular dynamics simulation with nuclear magnetic resonance (NMR) pseudocontact shift measurements. The 3site state is a result of simultaneous binding of three camphor molecules in three locations around P450cam: (a) in a well-established "catalytic" site near heme, (b) in a kink-separated "waiting" site along channel-1, and (c) in a previously reported "allosteric" site at E, F, G, and H helical junctions. These three spatially distinct binding modes in the 3site state mutually communicate with each other via homotropic allostery and act cooperatively to render P450cam functional. The 3site state shows a significantly superior fit with NMR pseudo contact shift (PCS) data with a Q-score of 0.045 than previously known bound states and consists of D251 free of salt-bridges with K178 and R186, rendering the enzyme functionally primed. To date, none of the reported cocomplex of P450cam with its redox partner putidaredoxin (pdx) has been able to match solution NMR data and controversial pdx-induced opening of P450cam's channel-1 remains a matter of recurrent discourse. In this regard, inclusion of pdx to the 3site state is able to perfectly fit the NMR PCS measurement with a Q-score of 0.08 and disfavors the pdx-induced opening of channel-1, reconciling previously unexplained remarkably fast hydroxylation kinetics with a koff of 10.2 s-1. Together, our findings hint that previous experimental observations may have inadvertently captured the 3site state as an in vitro solution state, instead of the catalytic state alone, and provided a distinct departure from the conventional understanding of cytochrome P450.

Redox partner recognition and selectivity of cytochrome P450lin (CYP111A1).

The strict requirement of cytochrome P450cam for its native ferredoxin redox partner, putidaredoxin (Pdx), is not exhibited by any other known cytochrome P450 (CYP) system and the molecular details of redox partner selectivity are still not completely understood. We therefore examined the selectivity of a related Pseudomonas cytochrome P450, P450lin, by testing its activity with non-native redox partners. We found that P450lin could utilize Arx, the native redox partner of CYP101D1, to enable turnover of its substrate, linalool, while Pdx showed limited activity. Arx exhibited a higher sequence similarity to P450lins native redox partner, linredoxin (Ldx) than Pdx, including several residues that are believed to be at the interface of the two proteins, based on the P450cam-Pdx complex structure. We therefore mutated Pdx to resemble Ldx and Arx and found that a double mutant, D38L/∆106, displayed higher activity than Arx. In addition, Pdx D38L/∆106 does not induce a low-spin shift in linalool bound P450lin but does destabilize the P450lin-oxycomplex. Together our results suggest that P450lin and its redox partners may form a similar interface to P450cam-Pdx, but the interactions that allow for productive turnover are different.

Hydroxylation Regiochemistry Is Robust to Active Site Mutations in Cytochrome P450cam (CYP101A1).

Cytochrome P450cam (CYP101A1) catalyzes the hydroxylation of d-camphor by molecular oxygen. The enzyme-catalyzed hydroxylation exhibits a high degree of regioselectivity and stereoselectivity, with a single major product, d-5-exo-hydroxycamphor, suggesting that the substrate is oriented to facilitate this specificity. In previous work, we used an elastic network model and perturbation response scanning to show that normal deformation modes of the enzyme structure are highly responsive not only to the presence of a substrate but also to the substrate orientation. This work examines the effects of mutations near the active site on substrate localization and orientation. The investigated mutations were designed to promote a change in substrate orientation and/or location that might give rise to different hydroxylation products, while maintaining the same carbon and oxygen atom balances as in the wild type (WT) enzyme. Computational experiments and parallel in vitro site-directed mutations of CYP101A1 were used to examine reaction products and enzyme activity. 1H-15N TROSY-HSQC correlation maps were used to compare the computational results with detectable perturbations in the enzyme structure and dynamics. We found that all of the mutant enzymes retained the same regio- and stereospecificity of hydroxylation as the WT enzyme, with varying degrees of efficiency, which suggests that large portions of the enzyme have been subjected to evolutionary pressure to arrive at the appropriate sequence-structure combination for efficient 5-exo hydroxylation of camphor.

Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution.

Steered molecular dynamics (SMD) simulation is a powerful method in computer-aided drug design as it can be used to access the relative binding affinity with high precision but with low computational cost. The success of SMD depends on the choice of the direction along which the ligand is pulled from the receptor-binding site. In most simulations, the unidirectional pathway was used, but in some cases, this choice resulted in the ligand colliding with the complex surface of the exit tunnel. To overcome this difficulty, several variants of SMD with multidirectional pulling have been proposed, but they are not completely devoid of disadvantages. Here, we have proposed to determine the direction of pulling with a simple scoring function that minimizes the receptor-ligand interaction, and an optimization algorithm called differential evolution is used for energy minimization. The effectiveness of our protocol was demonstrated by finding expulsion pathways of Huperzine A and camphor from the binding site of Torpedo California acetylcholinesterase and P450cam proteins, respectively, and comparing them with the previous results obtained using memetic sampling and random acceleration molecular dynamics. In addition, by applying this protocol to a set of ligands bound with LSD1 (lysine specific demethylase 1), we obtained a much higher correlation between the work of pulling force and experimental data on the inhibition constant IC50 compared to that obtained using the unidirectional approach based on minimal steric hindrance.

Updating the Paradigm: Redox Partner Binding and Conformational Dynamics in Cytochromes P450.

This Account summarizes recent findings centered on the role that redox partner binding, allostery, and conformational dynamics plays in cytochrome P450 proton coupled electron transfer. P450s are one of Nature's largest enzyme families and it is not uncommon to find a P450 wherever substrate oxidation is required in the formation of essential molecules critical to the life of the organism or in xenobiotic detoxification. P450s can operate on a remarkably large range of substrates from the very small to the very large, yet the overall P450 three-dimensional structure is conserved. Given this conservation of structure, it is generally assumed that the basic catalytic mechanism is conserved. In nearly all P450s, the O2 O-O bond must be cleaved heterolytically enabling one oxygen atom, the distal oxygen, to depart as water and leave behind a heme iron-linked O atom as the powerful oxidant that is used to activate the nearby substrate. For this process to proceed efficiently, externally supplied electrons and protons are required. Two protons must be added to the departing O atom while an electron is transferred from a redox partner that typically contains either a Fe2S2 or FMN redox center. The paradigm P450 used to unravel the details of these mechanisms has been the bacterial CYP101A1 or P450cam. P450cam is specific for its own Fe2S2 redox partner, putidaredoxin or Pdx, and it has long been postulated that Pdx plays an effector/allosteric role by possibly switching P450cam to an active conformation. Crystal structures, spectroscopic data, and direct binding experiments of the P450cam-Pdx complex provide some answers. Pdx shifts the conformation of P450cam to a more open state, a transition that is postulated to trigger the proton relay network required for O2 activation. An essential part of this proton relay network is a highly conserved Asp (sometimes Glu) that is known to be critical for activity in a number of P450s. How this Asp and proton delivery networks are connected to redox partner binding is quite simple. In the closed state, this Asp is tied down by salt bridges, but these salt bridges are ruptured when Pdx binds, leaving the Asp free to serve its role in proton transfer. An alternative hypothesis suggests that a specific proton relay network is not really necessary. In this scenario, the Asp plays a structural role in the open/close transition and merely opening the active site access channel is sufficient to enable solvent protons in for O2 protonation. Experiments designed to test these various hypotheses have revealed some surprises in both P450cam and other bacterial P450s. Molecular dynamics and crystallography show that P450cam can undergo rather significant conformational gymnastics that result in a large restructuring of the active site requiring multiple cis/trans proline isomerizations. It also has been found that X-ray driven substrate hydroxylation is a useful tool for better understanding the role that the essential Asp and surrounding residues play in catalysis. Here we summarize these recent results which provide a much more dynamic picture of P450 catalysis.

A Conserved Allosteric Site on Drug-Metabolizing CYPs: A Systematic Computational Assessment.

Cytochrome P450 enzymes (CYPs) are the largest group of enzymes involved in human drug metabolism. Ligand tunnels connect their active site buried at the core of the membrane-anchored protein to the surrounding solvent environment. Recently, evidence of a superficial allosteric site, here denoted as hotspot 1 (H1), involved in the regulation of ligand access in a soluble prokaryotic CYP emerged. Here, we applied multi-scale computational modeling techniques to study the conservation and functionality of this allosteric site in the nine most relevant mammalian CYPs responsible for approximately 70% of drug metabolism. In total, we systematically analyzed over 44 µs of trajectories from conventional MD, cosolvent MD, and metadynamics simulations. Our bioinformatic analysis and simulations with organic probe molecules revealed the site to be well conserved in the CYP2 family with the exception of CYP2E1. In the presence of a ligand bound to the H1 site, we could observe an enlargement of a ligand tunnel in several members of the CYP2 family. Further, we could detect the facilitation of ligand translocation by H1 interactions with statistical significance in CYP2C8 and CYP2D6, even though all other enzymes except for CYP2C19, CYP2E1, and CYP3A4 presented a similar trend. As the detailed comprehension of ligand access and egress phenomena remains one of the most relevant challenges in the field, this work contributes to its elucidation and ultimately helps in estimating the selectivity of metabolic transformations using computational techniques.

Partial Opening of Cytochrome P450cam (CYP101A1) Is Driven by Allostery and Putidaredoxin Binding.

Cytochrome P450cam (CYP101A1) catalyzes the regio- and stereo-specific 5-exo-hydroxylation of camphor via a multistep catalytic cycle that involves two-electron transfer steps, with an absolute requirement that the second electron be donated by the ferrodoxin, putidaredoxin (Pdx). Whether P450cam, once camphor has bound to the active site and the substrate entry channel has closed, opens up upon Pdx binding, during the second electron transfer step, or it remains closed is still a matter of debate. A potential allosteric site for camphor binding has been identified and postulated to play a role in the binding of Pdx. Here, we have revisited paramagnetic NMR spectroscopy data and determined a heterogeneous ensemble of structures that explains the data, provides a complete representation of the P450cam/Pdx complex in solution, and reconciles alternative hypotheses. The allosteric camphor binding site is always present, and the conformational changes induced by camphor binding to this site facilitates Pdx binding. We also determined that the state to which Pdx binds comprises an ensemble of structures that have features of both the open and closed state. These results demonstrate that there is a finely balanced interaction between allosteric camphor binding and the binding of Pdx at high camphor concentrations.

Controlling the Substrate Specificity of an Enzyme through Structural Flexibility by Varying the Salt-Bridge Density.

Many enzymes, particularly in one single family, with highly conserved structures and folds exhibit rather distinct substrate specificities. The underlying mechanism remains elusive, the resolution of which is of great importance for biochemistry, biophysics, and bioengineering. Here, we performed a neutron scattering experiment and molecular dynamics (MD) simulations on two structurally similar CYP450 proteins; CYP101 primarily catalyzes one type of ligands, then CYP2C9 can catalyze a large range of substrates. We demonstrated that it is the high density of salt bridges in CYP101 that reduces its structural flexibility, which controls the ligand access channel and the fluctuation of the catalytic pocket, thus restricting its selection on substrates. Moreover, we performed MD simulations on 146 different kinds of CYP450 proteins, spanning distinct biological categories including Fungi, Archaea, Bacteria, Protista, Animalia, and Plantae, and found the above mechanism generally valid. We demonstrated that, by fine changes of chemistry (salt-bridge density), the CYP450 superfamily can vary the structural flexibility of its member proteins among different biological categories, and thus differentiate their substrate specificities to meet the specific biological needs. As this mechanism is well-controllable and easy to be implemented, we expect it to be generally applicable in future enzymatic engineering to develop proteins of desired substrate specificities.

Covalent conjugation of single-walled carbon nanotube with CYP101 mutant for direct electrocatalysis.

Covalent linkage between the single-walled carbon nanotube (SWCNT) and CYP101 through a specific site of the enzyme can provide a novel method of designing efficient enzyme electrodes using this prototype cytochrome P450 enzyme. We have chemically modified the SWCNT with linker 4-carboxy phenyl maleimide (CPMI) containing maleimide functional groups. The enzyme was covalently attached on to the SWCNT through the maleimide group of the linker (CPMI) to the thiolate group of the surface exposed Cys 58 or Cys 136 of the CYP101 forming a covalently immobilized protein on the nanotube. Thin film of the modified SWCNT-CPMI-CYP101conjugate was made on a glassy carbon (GC) electrode. Direct electrochemistry of the substrate (camphor)-bound enzyme was studied using this immobilized enzyme electrode system and the redox potential was found to be -320mV vs Ag/AgCl (3 M KCl), which agrees with the redox potential of the substrate bound enzyme reported earlier. The electrochemically driven enzymatic mono-oxygenation of camphor by this immobilized enzyme electrode system was studied by measurement of the catalytic current at different concentrations of camphor. The catalytic current was found to increase with increasing concentration of camphor in presence of oxygen. The product formed during the catalysis was identified by mass-spectrometry as hydroxy-camphor.

Active Site Hydrogen Bonding Induced in Cytochrome P450cam by Effector Putidaredoxin.

Cytochrome P450s are diverse and powerful catalysts that can activate molecular oxygen to oxidize a wide variety of substrates. Catalysis relies on effective uptake of two electrons and two protons. For cytochrome P450cam, an archetypal member of the superfamily, the second electron must be supplied by the redox partner putidaredoxin (Pdx). Pdx also plays an effector role beyond electron transfer, but after decades the mechanism remains under investigation. We applied infrared spectroscopy to heme-ligated CN- to examine the influence of Pdx binding. The results indicate that Pdx induces the population of a conformation wherein the CN- ligand forms a strong hydrogen bond to a solvent water molecule, experimentally corroborating the formation of a proposed proton delivery network. Further, characterization of T252A P450cam implicates the side chain of Thr252 in regulating the population equilibrium of hydrogen-bonded states within the P450cam/Pdx complex, which could underlie its role in directing activated oxygen toward product formation and preventing reaction uncoupling through peroxide release.

Unexpected Differences between Two Closely Related Bacterial P450 Camphor Monooxygenases.

The bacterial cytochrome P450cam catalyzes the oxidation of camphor to 5-exo-hydroxycamphor as the first step in the oxidative assimilation of camphor as a carbon/energy source. CYP101D1 is another bacterial P450 that catalyzes the same reaction. A third P450 (P450tcu) has recently been discovered that has ≈86% sequence identity to P450cam as well as very similar enzymatic properties. P450tcu, however, exhibits three unusual features not found in P450cam. First, we observe product in at least two orientations in the X-ray structure that indicates that, unlike the case for P450cam, X-ray-generated reducing equivalents can drive substrate hydroxylation in crystallo. We postulate, on the basis of molecular dynamics simulations, that greater flexibility in P450tcu enables easier access of protons to the active site and, together with X-ray driven reduction, results in O2 activation and substrate hydroxylation. Second, the characteristic low-spin to high-spin transition when camphor binds occurs immediately with P450cam but is very slow in P450tcu. Third, isothermal titration calorimetry shows that in P450cam substrate binding is entropically driven with a ΔH of >0 while in P450tcu with a ΔH of <0 with a more modest change in -TΔS. These results indicate that despite nearly identical structures and enzymatic properties, these two P450s exhibit quite different properties most likely related to differences in conformational dynamics.

Proton Relay Network in the Bacterial P450s: CYP101A1 and CYP101D1.

Cytochrome P450s are among nature's most powerful catalysts. Their ability to activate molecular dioxygen to form high-valent ferryl intermediates (Compounds I and II) enables a wide array of chemistries ranging from simple epoxidations to more complicated C-H bond oxidations. Oxygen activation is achieved by reduction of the ferrous dioxygen complex, which requires the transfer of an electron from a redox partner and subsequent double protonation to yield a water molecule and a ferryl porphyrin π-cation radical (Compound I). Previous studies of the CYP101 family of cytochrome P450s demonstrated the importance of the conserved active site Asp25X residue in this protonation event, although its precise role is yet to be unraveled. To further explore the origin of protons in oxygen activation, we analyzed the effects of an Asp to Glu mutation at the 25X position in P450cam and in CYP101D1. This mutation inactivates P450cam but not CYP101D1. A series of mutagenic, crystallographic, kinetic, and molecular dynamics studies indicate that this mutation locks P450cam into a closed, inactive conformation. In CYP101D1, the D259E mutant changes the rate-limiting step to reduction of the P450-oxy complex, thus opening a window into the critical proton-coupled electron transfer step in P450 catalysis.

Proton relay network in P450cam formed upon docking of putidaredoxin.

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key O─O cleavage step remain partly elusive to date; effects observed in various enzyme mutants remain partly unexplained. We have carried out extended (to 1000 ns) MM-MD and follow-on quantum mechanics/molecular mechanics computations, both on the well-studied FeOO state and on Cpd(0) (compound 0). Our simulations include (all camphor-bound): (a) WT (wild type), FeOO state. (b) WT, Cpd(0). (c) Pdx (Putidaredoxin, redox partner of P450)-docked-WT, FeOO state. (d) Pdx-docked WT, Cpd(0). (e) Pdx-docked T252A mutant, Cpd(0). Among our key findings: (a) Effect of Pdx docking appears to go far beyond that indicated in prior studies: it leads to specific alterations in secondary structure that create the crucial proton relay network. (b) Specific proton relay networks we identify are: FeOO(H)⋯T252⋯nH 2 O⋯D251 in WT; FeOO(H)⋯nH 2 O⋯D251 in T252A mutant; both occur with Pdx docking. (c) Direct interaction of D251 with -FeOOH is, respectively, rare/frequent in WT/T252A mutant. (d) In WT, T252 is in the proton relay network. (e) Positioning of camphor appears significant: when camphor is part of H-bonding network, second protonation appears to be facilitated.

An Intermediate Conformational State of Cytochrome P450cam-CN in Complex with Putidaredoxin.

Cytochrome P450cam is an archetypal example of the vast family of heme monooxygenases and serves as a model for an enzyme that is highly specific for both its substrate and reductase. During catalysis, it undergoes significant conformational changes of the F and G helices upon binding its substrate and redox partner, putidaredoxin (Pdx). Recent studies have shown that Pdx binding to the closed camphor-bound form of ferric P450cam results in its conversion to a fully open state. However, during catalytic turnover, it remains unclear whether this same conformational change also occurs or whether it is coupled to the formation of the critical compound I intermediate. Here, we have examined P450cam bound simultaneously by camphor, CN-, and Pdx as a mimic of the catalytically competent ferrous oxy-P450cam-Pdx state. The combined use of double electron-electron resonance and molecular dynamics showed direct observation of intermediate conformational states of the enzyme upon CN- and subsequent Pdx binding. This state is coupled to the movement of the I helix and residues at the active site, including Arg-186, Asp-251, and Thr-252. These movements enable occupation of a water molecule that has been implicated in proton delivery and peroxy bond cleavage to give compound I. These findings provide a detailed understanding of how the Pdx-induced conformational change may sequentially promote compound I formation followed by product release, while retaining stereoselective hydroxylation of the substrate of this highly specific monooxygenase.

Cytochrome P450-The Wonderful Nanomachine Revealed through Dynamic Simulations of the Catalytic Cycle.

This Account addresses the catalytic cycle of the enzyme cytochrome P450 (CYP450) as a prototypical biological machine with automatic features. CYP450 is a nanomachine that uses dioxygen and two reducing and two proton equivalents to oxidize a plethora of molecules (so-called substrates) as a means of supplying bio-organisms with essential molecules (e.g., brain neurotransmitters, sex hormones, etc.) and protecting biosystems against poisoning. An enticing property of CYP450s is that entrance of an oxidizable substrate into the active site initiates a series of events that constitute the catalytic cycle, which functions "automatically" in a regulated sequence of events culminating in the production of the oxidized substrates (e.g., hydroxylated, epoxidized, etc.), oftentimes with remarkable stereo- and regioselectivities. It is timely to demonstrate how theory uses molecular dynamics (MD) simulations and quantum-mechanical/molecular-mechanical (QM/MM) calculations to complement experiments and elucidate the choreography by which the protein regulates the catalytic cycle. CYP450 is a heme enzyme that contains a ferric ion (FeIII) coordinated by a porphyrin ligand, a water molecule, and a cysteinate ligand that is provided by a strategic residue of the encapsulating protein. While many of the individual steps are sufficiently well-understood, we shall provide here an overview of the factors that cause all of the steps to be sequentially coordinated. To this end, we use examples from three different CYP450 enzymes: the bacterial ones CYP450BM3 and CYP450CAM and the mammalian enzyme CYP4503A4. The treatment is limited to the catalytic cycle, as aspects of two-state reactivity were reviewed previously (e.g., Shaik , S. ; et al. Chem. Rev. 2005 , 105 , 2279 ). What are the principles that govern the seeming automatic feature? For example, how do substrate entrance and binding gate the enzyme? How does the reductase attachment to the enzyme affect the next steps? What triggers the attachment of the reductase? How does the electron transfer (ET) that converts FeIII to FeII occur? Is the ET coordinated with the entrance of O2 into the active site? What is the mechanism of the latter step? Since the entrance of the substrate expels the water molecules from the active site, how do water molecules re-enter to form a proton channel, which is necessary for creating the ultimate oxidant Compound I? How do mutations that disrupt the water channel nevertheless create a competent oxidant? By what means does the enzyme produce regio- and stereoselective oxidation products? What triggers the departure of the oxidized product, and how does the exit occur in a manner that generates the resting state ready for the next cycle? This Account shows that the entrance of the substrate triggers all of the ensuing events.

Stereo-preference of camphor for H-bonding with phenol, methanol and chloroform: A combined matrix isolation IR spectroscopic and quantum chemical investigation.

Camphor is known to be held in the substrate pocket of cytochrome P450cam enzyme via H-bond with a tyrosine residue of the enzyme in a unique orientation. This structural exclusivity results in regio- and stereo-specific hydroxylation of camphor by the enzyme. We have carried out a combined IR spectroscopic and quantum chemical investigation to shed light on the factors influencing the conformational exclusivity of 1R-(+)-camphor in the substrate pocket of Cytochrome P450cam, and to determine whether the selectivity is an inherent property of the substrate itself, or is imposed by the enzyme. For this purpose, complexes of camphor have been studied with three H-bond donors namely phenol, methanol and chloroform. Each of the three donors was found to form stable complexes with two distinct conformers; the one mimicking the conformation in enzyme substrate pocket was found to be more stable of the two, for all three donors. Experimentally, both conformers of the H-bonded complexes were identified separately for phenol and methanol in an argon matrix at 8 K, but not for chloroform due to very small energy barrier for interconversion of the two conformers. In room temperature solution phase spectra of camphor with all three donors, the differences in spectral attributes between the two isomeric H-bonded complexes were lost due to thermal motions.

Mapping the Substrate Recognition Pathway in Cytochrome P450.

Cytochrome P450s are ubiquitous metalloenzymes involved in the metabolism and detoxification of foreign components via catalysis of the hydroxylation reactions of a vast array of organic substrates. However, the mechanism underlying the pharmaceutically critical process of substrate access to the catalytic center of cytochrome P450 is a long-standing puzzle, further complicated by the crystallographic evidence of a closed catalytic center in both substrate-free and substrate-bound cytochrome P450. Here, we address a crucial question whether the conformational heterogeneity prevalent in cytochrome P450 translates to heterogeneous pathways for substrate access to the catalytic center of these metalloenzymes. By atomistically capturing the full process of spontaneous substrate association from bulk solvent to the occluded catalytic center of an archetypal system P450cam in multi-microsecond-long continuous unbiased molecular dynamics simulations, we here demonstrate that the substrate recognition in P450cam always occurs through a single well-defined dominant pathway. The simulated final bound pose resulting from these unguided simulations is in striking resemblance with the crystallographic bound pose. Each individual binding trajectory reveals that the substrate, initially placed at random locations in bulk solvent, spontaneously lands on a single key channel on the protein-surface of P450cam and resides there for an uncharacteristically long period, before correctly identifying the occluded target-binding cavity. Surprisingly, the passage of substrate to the closed catalytic center is not accompanied by any large-scale opening in protein. Rather, the unbiased simulated trajectories (∼57 µs) and underlying Markov state model, in combination with free-energy analysis, unequivocally show that the substrate recognition process in P450cam needs a substrate-induced side-chain displacement coupled with a complex array of dynamical interconversions of multiple metastable substrate conformations. Further, the work reconciles multiple precedent experimental and theoretical observations on P450cam and establishes a comprehensive view of substrate-recognition in cytochrome P450 that only occurs via substrate-induced structural rearrangements.

Substrate-Dependent Allosteric Regulation in Cytochrome P450cam (CYP101A1).

Various biophysical methods have provided evidence of a second substrate binding site in the well-studied cytochrome P450cam, although the location and biological relevance of this site has remained elusive. A related question is how substrate and product binding and egress occurs. While many active site access channels have been hypothesized, only one, channel 1, has been experimentally validated. In this study, molecular dynamics simulations reveal an allosteric site related to substrate binding and product egress. The remote allosteric site opens channel 1 and primes the formation of a new channel that is roughly perpendicular to channel 1. Substrate entry to the active site via channel 1 as well as substrate/product egress via channel 2 is observed after binding of a second molecule of substrate to the allosteric site, indicating cooperativity between these two sites. These results are consistent with and bring together many early and recent experimental results to reveal a dynamic interplay between a weak allosteric site and substrate binding to the active site that controls P450cam activity.

Direct Experimental Characterization of Contributions from Self-Motion of Hydrogen and from Interatomic Motion of Heavy Atoms to Protein Anharmonicity.

One fundamental challenge in biophysics is to understand the connection between protein dynamics and its function. Part of the difficulty arises from the fact that proteins often present local atomic motions and collective dynamics on the same time scales, and challenge the experimental identification and quantification of different dynamic modes. Here, by taking lyophilized proteins as the example, we combined deuteration technique and neutron scattering to separate and characterize the self-motion of hydrogen and the collective interatomic motion of heavy atoms (C, O, N) in proteins on the pico-to-nanosecond time scales. We found that hydrogen atoms present an instrument-resolution-dependent onset for anharmonic motions, which can be ascribed to the thermal activation of local side-group motions. However, the protein heavy atoms exhibit an instrument-resolution-independent anharmonicity around 200 K, which results from unfreezing of the relaxation of the protein structures on the laboratory equilibrium time (100-1000 s), softening of the entire bio-macromolecules.

Effect of redox partner binding on CYP101D1 conformational dynamics.

We have compared the thermodynamics of substrate and redox partner binding of P450cam to its close homologue, CYP101D1, using isothermal titration calorimetry (ITC). CYP101D1 binds camphor about 10-fold more weakly than P450cam which is consistent with the inability of camphor to cause a complete low- to high-spin shift in CYP101D1. Even so molecular dynamics simulations show that camphor is very stable in the CYP101D1 active site similar to P450cam. ITC data on the binding of the CYP101D1 ferredoxin redox partner (abbreviated Arx) shows that the substrate-bound closed state of CYP101D1 binds Arx more tightly than the substrate-free open form. This is just the opposite to P450cam where Pdx (ferredoxin redox partner of P450cam) favors binding to the P450cam open state. In addition, CYP101D1-Arx binding has a large negative ΔS while the P450cam-Pdx has a much smaller ΔS indicating that interactions at the docking interface are different. The most obvious difference is that PDXD38 which forms an important ion pair with P450camR112 at the center of the interface is ArxL39 in Arx. This suggests that Arx may adopt a different orientation than Pdx in order to optimize nonpolar interactions with ArxL39.